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ABSTRACT 

We have used the rocket-borne Wide-Field Imaging Survey Polarimeter (WISP) to image a 1?5 x 4?8 
area of the western side of the Large Magellanic Cloud (LMC) at a wavelength of A = 2150 A and a 
resolution of 1' x 1.'5. These are the first wide- field ultraviolet polarimetric images in astronomy. We 
find the UV background light of the LMC to be linearly polarized at levels ranging from our sensitivity 
limit of 4% to as high as «40%. In general, the polarization in a pixel increases as the flux decreases; 
the weighted mean value of polarization across the WISP field is 12.6% ±2.3%. The LMC's diffuse UV 
background, in uncrowded areas, rises from a minimum of 5.6 ±3.1 x 1CU 8 erg s _1 cm~ 2 A -1 Sr _1 
(23.6 ±0.5 mag arcsec" 2 ) to 9.3 ±1.1 x 10~ 8 erg s" 1 cm" 2 A" 1 Sr" 1 (23.1 ±0.2 mag arcsec -2 ) in 
regions near the bright associations. We use our polarization maps to investigate the geometry of the 
interstellar medium in the LMC, and to search for evidence of a significant contribution of scattered 
light from OB associations to the diffuse galactic light of the LMC. Through a statistical analysis of our 
polarization map, we identify 9 regions of intense UV emission which may be giving rise to scattering 
halos in our image. We find that starlight from the N 11 complex and the LH 15 association are the 
strongest contributors to the scattered light component of the LMC's diffuse galactic light. This region 
of the northwestern LMC can be thought of as a kiloparsec-scale reflection nebula in which OB stars 
illuminate distant dust grains, which scatter the light into our sightline. In contrast, the polarization 
map does not support the scattering of light from the large B2 complex in the southern WISP field; this 
effect may be astrophysical, or it may be the result of bias in our analysis. 

Subject headings: galaxies: individual (LMC) — polarization — ultraviolet emission — stars: 
early-type — ISM: structure 
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1. INTRODUCTION 



to its depletion in carbo n relative to the Milky Way; Pe 



rpi r i j j i- i r i i i l j • 1992, Fitzpatrick (1986)) should compensate somewhat. 
1 no presence of extended optical light halos around spi r ^ — . ,\ (V Q f . /wTcm • 



ral galaxies raises the question: is this light due to an 
extended stellar distribution, or is it the re sult of scat- 
tered starlight from more central regions? |Jura (1980)| 
pointed out that if the extended clouds of neutral hydro- 
gen around late- type galaxies are associated with dust, the 
clouds should appear as faint reflection nebulae. In their 
investigation of the diffus e ultraviolet light of the Galaxy, 
Murthy & Henry (1995) found a significant contribution 
due to the scattering of starlight from OB stars by inter- 
stellar dust. Photometric studies of spiral galaxies have 
provided very suggestive evidence for UV scattering ha- 
los far from reg ions of star-formation (Stecher et al. 1982 
Ncff et al. 1994). Signatures of dust-scattered starlight are 



spatial correlations with H I, bluer colors than the general 
radiation field, and the presence of linear polarization. 
The LMC, due to its proximity, vigorous star-forming 



activity, and low level of obscuration by Galactic du3t, is 



an ideal location in which to investigate the connection 

between starlight and diffuse light in a spiral galaxy. The 
LMC is known to possess a higher gas -to-dust ratio than 
the Galaxy (e.g., Clayton et al. 1996); while this would 



tend to lessen the contribution of scattering to its diffuse 
galactic light, the relatively high albedo of LMC dust (due 



The Wide-Field Imaging Survey Polarimeter (WISP) in- 
strument is designed to obtain polarimetric images of as- 
tronomical objects in the ultraviolet in order to investigate 
the composition and distribution of dust in astrophysical 
environments. The ultraviolet is chosen because the emis- 
sion from late-type stars, which are smoothly distributed 
and numerous, is suppressed. At optical wavelengths, late- 
type stars contribute strongly to the general radiation field 
of a galaxy and may mask the presence of scattered light 
with their own radiation. In the ultraviolet, OB stars in 
spatia lly compact associati ons dominate the radiation field 



(e.g., Parker et al. 199S for the LMC); far from these 



associations, scattered light at these wavelengths is ex- 
pected to contribute strongly to the diffuse radiation field. 
In addition to the low stellar background, the UV is an 
ideal regime in which to make polarimet ric investigations 



becau se of high polarimetric efficiency (Nordsieck et al 



1993) 



If scattering is a significant contributor to the diffuse ra- 
diation field, then its surface brightness and polarimetric 
properties can be used to derive information regarding the 
distribution of dust in the interstellar medium, as well as 
to place constraints on the optical properties of the dust 
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grains. 

WISP has been flown three times to date; its targets 
have be en the reflection nebulosity in the Plei a des star 
cluster ( |Gibson, Holdaway, fc Nordsicck 1995 ; |Gibson| 
19970 Comet Hale-Bopp ( Harris et al. 19971 ), and tnc 
LMC. 

In this paper, we present polarimetric images of a field 
along the western side of the LMC. These are the first 
wide-field ultraviolet polarimetric images ever obtained in 
astronomy; ground-based imaging polarimetry in the rest 
frame UV of radio galaxie s and QSOs at redsh i fts z ~ 1 



has b een analyzed by, e .g., |Cimatti et al. (1994) 
Aligheri et al. (1993), and references therein. Capetti ei 



|di~STr 



polarization Stokes parameters: "Q + " is obtained with 
zero waveplate pressure; "Q~", "U~", and "U + " with i 
wave retardation, rotated 45°, 67?5, and 22?5 relative to 
the Brewster mirror axis, respectively. Because the mod- 
ulator is the first clement in the WISP optical path, in- 
strumental polarization from subsequent optical elements 
is eliminated. The position angle zero-point was fixed by 
optical alignment and mechanical tolerances to better than 
±1°. 

Polarimetric calibration was performed in August, 1995 
using the 40 cm vacuum collimator at the Space Astron- 
omy Lab, University of Wisconsin. 100% linearly polarized 
continuum ultraviolet light was created using a D2 hollow 
cathode lamp illuminating the focal plane aperture of the 
collimator, immediately followed by an MgF2 Brewster- 
angle polarizer. The polarimetric efficiency of the analyzer 
is estimated to be better than 98%, based on crossing two 
analyzers in visible light, and on analysis of the effect of 
small changes in the Brewster angle in the vacuum ultra- 
violet. Collimated light from this device illuminated the 
center and extrema of the WISP field by maneuvering the 
payload on a remotely controlled tilt table. The efficiency 
of the WISP analyzer was evaluated by rotating the cali- 
bration polarizer with the waveplate off; it is 94% in the 
2150 A filter. The optimal stress for halfwave retardation 
was evaluated by fixing the calibration polarizer angle and 
WISP waveplate angle, while varying the waveplate stress. 
The waveplate efficiency was then evaluated by modeling 
the polarimetric modulation with waveplate angle at the 
optimal stress, and found to be 79% in the 2150 A filter. 
The net polarimetric efficiency in this filter is therefore 
«74%; a figure of 70% was used during reductions. 

3. DATA 

3.1. Flight and Observations 

WISP was flown from Woomera, Australia on 17 
November, 1995; the flight was a complete success. For the 
LMC flight the 15" CCD pixels were binned to 4 columns 
by 6 rows, resulting in an image scale of 1' x 1.'5 per pixel. 
The 2150 A filter was used to take a polarimetric sequence 
of 4 80-second exposures, in the order Q~, Q + , U~, U + . 
The science exposures began at an altitude of 230 km, con- 
tinued through the apogee of 328 km, and ended at 150 
km. The planned pointing placed the long axis of the CCD 
9° from north-south, with the field centered at a = 4 h 58 m 
45 s , S = -67° 55' (B1950.0). The pointing drifted by less 
than 5' between the first and last exposures. 

3.2. Reduction 

We reduced the images within IRAF 1 , removing the bias 
and dark count in the usual way. During the flight, a se- 
quence of 8 bias frames and 3 40-second dark frames were 
taken for calibration. Additionally, a sequence of 12 dark 
exposures taken in the laboratory were combined to pro- 
duce a superdark frame. The lab darks confirmed the lin- 
earity of the CCD, and enabled us to better characterize 
the small-scale variations in pixel response when subtract- 
ing the dark frames. 

The rocket had begun its descent by the time of the last 
exposure, and a significant additional sky background was 

X IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



al. (1)95) used the Hubble Space Telescope's Faint Object 
Camera to obtain a UV polarimetric image of the galaxy 
NGC 1068 with an 11" field of view. 

We describe the WISP instrument in §||; in §|| we pro- 
vide the flight information and observation sequence, and 
describe the procedures we have adopted for reduction and 
calibration of the images. §|| gives the results of our sur- 
face photometry of the diffuse galactic light of the LMC. In 
§H we present our polarization map and our analysis which 
reveals the scattering halos around some of the prominent 
concentrations of young stars. In §|| we analyze the mor- 
phology of the UV image and polarization map and how it 
fits in with what is known about the structure and content 
of the ISM in this region. We discuss in §[7] the origin of the 
diffuse UV radiation field in the LMC and the scattering 
geometry within the LMC disk; our results are summa- 
rized in §||. 

2. THE INSTRUMENT 

The Wide-field Imaging Survey Polarimeter (WISP) is 
a suborbital rocket payload based around a 20 cm, F/1.8, 
off-axis, all-reflective Schmidt telescope. Nordsieck et al. 
(1993) provides a detailed description of the payload op- 
tics. The Schmidt corrector mirror removes the spherical 
aberration of the primary, providing a 1?5 x 4? 8 degree 
field of view with resolution well-suited to the 15 arcsec- 
ond pixel scale of the CCD detector. The WISP detector 
for this flight was a thick, front-illuminated Reticon 400 x 
1200 pixel CCD, phosphor overcoated to provide UV sen- 
sitivity. The quantum efficiency of the detector in the vac- 
uum ultraviolet is approximately 15% and the readnoise 
is 10 e~. The CCD was cooled to —65 °C by a thermo- 
electric cooler which dissipates its heat into a copper block 
which is precooled to —40 °C on the ground. WISP has 
two broadband filters, at 1640 ±250 A and 2150 ±300 A, 
each comprised of aluminum/MgF2 thin films on MgF2 
and silica substrates. 

The polarimetric analyzer is a large flat mirror mounted 
between the Schmidt corrector and primary, coated with 
Zr02 (a high- index dielectric), and illuminated at its 
Brewster angle. The polarimetric modulator is a rotatable 
20-cm halfwave plate mounted in the telescope aperture. 
The waveplate is constructed from four panes of oriented 
CaF2 crystal. The crystal is made to be birefringent by 
applying stress to one side using a programmable actuator. 
The entire waveplate assembly is then rotated to provide 
polarimetric modulation. Four images provide the linear 
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present in the U + image. We corrected for this by requir- 
ing that the sum of the Q images be equal to the sum 
of the U images, subtratcting the difference in sky back- 
ground from the U + image. The correction amounted to 
~16 counts, some 7.7% of the total sky background. 

We removed cosmic ray contamination in a two-step 
process. First, we used the IRAF task cosmicrays to 
find and remove pixels which deviated from the median 
of their neighbors by greater than a factor of ten. These 
regions were replaced by the mean of their neighbors. The 
remaining cosmic rays were found during a trial combina- 
tion of the Stokes images, which rejected pxiels that were 
more than 10a deviant from the mean across all four im- 
ages. Deviant pixels were masked off and excluded from 
further analysis. 10 bright, compact areas of emission, 
well-distributed across the field of view, were used to align 
and trim the images. The reduced images were combined 
to create intensity, Stokes Q, and U images: 



I = 



• Q- + u+ + in 

Q = Q + - Q~ 
U = U+ - IT 



(1) 
(2) 
(3) 



The intensity image was compared a gainst the catalog 
of iLucke fc Hodge (1970)1, the atlas of [Hodge fc Wright 



(19671, the UV images o f pmith et at (1987f (hereinafter 
SCH), and the catalog of fjanduleak 197C in order to iden- 
tify emission sources. We asso ciated 21 bright regions 



with the OB associations from |Lucke fc Hodge (1970) 
Other features include young open clusters, and in some 
uncrow ded regions we d etect individual, early-type super- 
giants ( Sandulcak 1970] ). We placed the image on a world 
coordinate system using some of these identifications and 



foundithc long axis of the image to be aligned 8° east of 
north The grayscale mlonsily image is overplollod on an 
optical image of the LMC in Figure |l} A schematic figure 
depicting some of the emission sources is shown in Figure 



3.3. Photometric Calibration 

SCH presented phot ometry at 1500 A and 1900 A for 
each association in the [Lucke fc Hodge (1970)| catalog. We 
used their results to obtain a photometric calibration of 
our data, using the equation 



in 



21 



mig - (mig - m 2 i)o 



Eiq—o 



(4) 



where mig refers to the measured 1900 A magnitude from 
SCH, (wig — m2i)o refers to the intrinsic 1900—2150 A 
color of an association, and Eig_2i denotes the amount of 



interstellar extinction distinguish the two filter bandpasses 
from each other. 

The list of calibrating associations is given in Table 1, 
which includes the magnitude at 1900 A and 1500-1900 A 
color from SCH, the aper ture used by S CH, and the E#_y 



for each association from Lucke (1974) 



We dereddened the photometry of SCH, using the value 
of Eb~v from Table 1 and assuming that the maximum 
contribution of the Galactic foreground was 0.07 mag. 
Lucke (1974) does not give E#_y for the associations LH 



10, 11, 13, or 14 and so we have estimated the reddening to 



be 0.07. We determined Ai 50 o/E B 
conv olving the extinctio n curves of 



and AiWEb^v by 



Cardelli et al. (1989) 



and Fitzpatrick (1986), assuming Ry = 3.1, with the 



source spectrum of a T e g = 30,000 K star (Kurucz 1991) 
and the filter transmission curves published by SCH. The 
resulting reddening corrections were: 



"115,0 = m 15 



8.82 Eg^y 



8.23 E A B m v 



OS rpLMC o 01 J?MW 
.ZSti B _ v — 8.21 Hi B -V- 



(5) 
(6) 



We then used the calculated 



mi5_ 



color, along 

with stellar atmospheres from [Kurucz 1991| and the filter 
transmission curves of SCH and the WISP 2150 A filter, 
to predict the (mig — m2i)o color of each association. For 
OB stars, we adopted a linear relation between the two 
colors: 



(mig - 17121)0 = 0.57 (mis - m w ) . 



(7) 



This relation becomes invalid for sources cooler than T e ff 
= 12,000 K. The intrinsic color of the sources suffers ad- 
ditional reddening due to the well-known extin ction bump 
at 2175 A. Using the LMC extinction curve of Fitzpatrick 



(1986) and the Milky Way curve of Cardelli et al. (1989 



we calculate 



ElQ- 



19-21 



-En- 



B-V, 



(8) 

note that because of the 2175 A bump, (mig — m2i) be- 
comes bluer with increasing extinction 3 

Substituting Equations |7] and ^ into Eq. g, we obtain 



syn 



TOig - 0.57 (77745 - mig) + E B - 



(9) 



selective extinction between the two wavelengths, m^f" 
refers to the synthetic 2150 A magnitude 2 of each associa- 
tion, calculated from the quantities on the right-hand side 
of Eq. g. 

For a flat spectrum source with zero reddening, F\ = 
constant, and Eq. g becomes m^f™ = mig. The two correc- 
tive terms arise because the complex spectral energy distri- 
bution of the associations and the wavelength-dependent 

2 We adopt throughout the definition m\ = —2.5 log (F^) — 21.1. 

3 The adoption of a mean UV extinction curve is fraught with peril, as well-documented by Fitzpatrick (1998) 
curve from sightline to sightline may contribute an additional 20% systematic uncertainty in photometry. 



The results for each association are given in column 7 of 
Table 1. 

We photometered the associations in the WISP image, 
attempting to match the apertures to those used by SCH. 
SCH do not give the precise locations of their apertures, 
but they provide the area of each region; we matched 
the areas using rectangular apertures in the IRAF task 
polyphot. The instrumental magnitudes are given in col- 
umn 6 of Table 1. 

We rejected the association LH 20 because it fell along 
the edge of the field; there was some uncertainty in the 
position of LH 14. We calibrated our photometry using a 
least-squares fit to m 2 i ; allowing for a nonlinearity as well 



as zeropoint offset, we found 



syn 



1.05 ±0.10) 



l 21 



(9.33 ±0.24), with an rms scatter of 0.26 mag. However, 



variation in the extinction 
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Fig. 2. — Prominent features in the WISP field, after Lucke & Hodge (1970). Numbered ellipses represent the positions of OB associations 
from Lucke & Hodge; the dashed ellipses N 11 and B2 refer to the enclosed OB associations and emission nebulae. Massive open clusters 
are denoted by pentagons and labelled with their NGC number. Sources of scattered light from Table 2 are marked with stars; large stars 
indicate stronger detections. 
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this fit was only marginally better than one in which we 
enforced linearity, and which we decided to adopt: 



as well as the uncertainties in photometric calibration. 
Due to resolution effects, Figure contains more extreme 



syn 

hi — 



values, both high and low, than the Maucherat-Joubcrt 



'21 



.53 ±0.28), 



(10) 



et al. (1980) map; the effect was lessened, although not 



with an rms scatter of 0.29 mag. For comparison to m^", 
our calibrated photometry for the 20 measured associa- 
tions is given in column 8 of Table 1. LH 14 was excluded 
from the fit because of the uncertainty in its identification; 
it is f aint, omall, and oubjoct to uncertain oky [subtraction 



eliminated, by smoothing the WISP image to 1° resolu 
tion. The difference is greatest for the degree-scale arcs 
seen in the WISP image north of w —67° 30': they are 
blended out by the 3° resolution in th e earlier work. The 
arc containing N 11 is visible in the Maucherat-Joubert 



due its close proximity to the extremely bright associations 
LH 9 and 10. 

3.4. Background Subtraction 



UVi surface brightness measurements arc affe cted by 



et al. (1980)| map as a small perturbation to the 11 x 
10~ 8 erg s _1 cm~ 2 A -1 Sr _1 contour line. In the south- 
ern WISP field, the light is more smoothly distributed, and 
the 16 x 10~ 8 erg s _1 cm" 2 A" 1 Sr _1 conto ur in Fig 
matches the location of the similar isophote in Maucherat 



three ma jor Dourc oo of background contamination (Murthy 



Joubert et al. (1980) 



fc Henry (1995) ): zodiacal light, Galactic background 
light, and the extragalactic background. In order to mea- 
sure the diffuse galactic light of t he LMC, we must re- 
move thes e background sky levels. Maucherat-Joubcrt e\ 
al. (1980) surveyed the entire sky at 2200 A with 3° rcsolu- 
tion, and measured a mean background sky level (Galactic 
plus extragalactic) of 0.9 ±0.2 x 10~ 8 erg s _1 cm~ 2 A -1 
Sr _1 for the Galactic latitude of the western LMC (—35° 
< b < -38°). They reported a zodiacal light level at 2200 



A of 0.10 ±0.02 x 10- 



erg s 1 cm 2 



for the 



LMC's ecliptic latitude. We therefore adopt a total back- 
ground level of 1.0 ±0.2 x 10~ 8 erg s^ 1 cm" 2 A" 1 Sr" 1 
(25.4 ±0.3 mag arcsec" 2 ) across the WISP field. This 
value has been subtracted from each of our surface bright- 
ness measures of the LMC's diffuse galactic light. 

4. SURFACE PHOTOMETRY 

A contour plot of the LMC's surface brightness at 2150 
A is shown in Figure at a resolution of 0?1. The dif- 
fuse UV background, which we attempted to measure in 
uncrowded areas, rises from a minimum of 5.6 ±3.1 x 

Sr~ x (23.6 ±0.5 mag arcsec -2 ) 
• ..,.„ o-i A-i Sr -i (23.0 ±0.2 



10~ 8 erg s -1 cm -2 A" 1 



to 9.3 ±1.1 x 10~ 8 erg s -1 cur 2 A 



mag arcsec ) in regions near the bright associations. The 
lowest areas occur where the stellar density is low, at the 
northern edge of the image, along the western edge of the 
image near —68°, and in the far southwestern corner. Un- 
crowded regions of high surface brightness include the re- 
gion between N 11 and LH 19, and the area north of NGC 
1755. Regions above the 11 x 10~ 8 erg s _1 cm~ 2 A -1 
Sr -1 contour line are highly resolved in the SCH images 
into sheets of stars; Fig. || shows that these high surface 
brightness regions are clustered around OB associations. 

The eastern edge of the image, at roughly —69°, in- 
cludes the extreme western end of the LMC bar. The bar 
is prominent in optical images of the galaxy, less so in UV 
images, e.g., SCH. The stellar density is far higher here 
than in the northern half of the WISP field; despite the 
suppression of late-type stars by the UV filter, the densely 
populated bar contributes to the high UV surface bright- 
ness. As will be discussed in this will adversely affect 
our ability to detect the effects of du st scattering. 

Maucherat-Joubert et al. (lyHTT) 1 presented a surface 
brightness map of the LMC at 2200 A. Agreement between 
the WISP sur face photometry and the Maucherat-Joubert 
et al. (1980) results is surprisingly good, considering the 



5. ULTRAVIOLET POLARIZATION MAP 

We mapped the degree of linear polarization, P, and 
position angle, 9 for the WISP field using the I, Q, and 
U images. Because the direct light of stars contributes 
highly to the flux but is largely unpolarized, their pres- 
ence tends to bias the estimates of the polarization of the 
diffuse galactic light to low values. To remove them from 
the polarization map, we cleaned the images of pixels that 
deviated by more than 5<r above the local sky value in 
all three of the images. The star mask was applied to the 
individual images and masked pixels were subsequently ig- 
nored. 

In order to increase the signal-to-noise, we binned the 
image into 36 arcmin 2 pixels (4x6 pixels). For each pixel 
we calculated the degree of linear polarization: 



/ (Q 2 + U 2 )\ ' 5 



V i 2 



(11) 



and the position angle, measured clockwise from the pos- 
itive y-axis (roughly north): 



6 = 0.5 arctan 



U 



(12) 



disparity in angular resolution and detector contstruction, 



We assumed Poisson statistics for the calculation of up 
and o~q. Because P is the sum in quadrature of Q and U, 
the expectation value of P is biased such that for constant 
Q and U, the most probable measured value of P increases 
with increasing <tq and o~\j. Because our data suffer from 
poor signal-to- noise levels, we hav e removed this bias using 
the formula of Nordsieck (1976). By this formula, a 2.5er 
detection is biased so that P bs/Ptrue « 1.1. We reject as 
non-detections pixels in which P < 2.5<7p. This restriction 
set our detection limit at 4% polarization. 

The final polarization map is shown in Figure |]a. Po- 
larizations as high as 40% are observed in some pixels, 
although their absolute flux level is small. Most of the 
values lie between 5% and 20%, with an error-weighted 
mean P = 12.6% ±2.3%. If we relax the restriction on 
detections below 2.5cr, the mean P = 4.4% ±7.8%. 

The maximum predicted polarization from a commonly 
used model of the interstellar dust mixture (Mathis, 
Rumpl, & Nordsieck 1977) isp max = 0.3-0.4 (White 1979). 
We can estimate the fractional contribution of scattered 
light to the total UV surface brightness in the WISP field 




Fig. 3. — Surface brightness plot of the WISP field; the image has been smoothed to 0?1 resolution. North is up and East is to the left. 
Contour lines are drawn at 6, 11, 16, 21, 26 X 1CT 8 erg s" 1 cm" 2 A" 1 Sr" 1 (23.5, 22.9, 22.5, 22.2, 21.9 mag arcsec^ 2 ). 
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Fig. 4. — a) Ultraviolet polarization map of the WISP field, binned to a resolution of 6 ' per pixel. North is up and East is to the left. Our 
map becomes incomplete for polarizations below 10%, and our sensitivity limit is 4%. Hints of centro-symmetric scattering patterns can be 
seen, b) A smoothed contour map of x 2 , where the x 2 value of a pixel represents the likelihood that the polarized light in neighboring pixels 
consists of singly-scattered photons originating in the pixel of interest. Contours arc plotted at \ 2 = 1>2,3,4,5. 
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by comparing our mean P to p m ax ■ This comparison sug- 
gests that 10-15% of the UV surface brightness is due to 
polarized scattered light. However, polarization as high 
as p m ax occurs only for scattering at a 90° angle, and de- 
creases for smaller and larger scattering angles. Because 
of the forward-throwing nature of dust grains in the ultra- 
violet, the mean scattering angle is expected to be much 
smaller than 90° for our data. The polarization of a scat- 
tered photon decreases asymptotically to zero as the scat- 
tering angle approaches zero (e.g., Scarrott et al. 1990). 
We therefore predict that much more than 10-15% of the 
measured UV surface brightness is contributed by scat- 
tered photons. 

For the case of scattering from a point source, the po- 
larization vectors form a characteristic centro-symmetric 
pattern around the illuminating source. The typical er- 
ror in our position angles is ±25°; this noise distorts 
the centro-symmetric signature of scattering halos around 
bright sources, if such features are present. Although our 
map is noisy, the high mean level of P is indicative of a 
significant scattered light component to the general UV 
surface brightness. 

To identify the sources most likely responsible for il- 
luminating the interstellar dust, we performed a statis- 



tical analysis designed to evaluate the degree of centro- 
symmetry m the polarization map. We stepped through 
the polarization map, pixel-by-pixel, and for each pixel, 
we examined the distribution of polarization vectors in 
the encompassing 7x7 pixel (42' x 42') region. For each 
region, we computed a x 2 value based on the distribution 
of (6 O0S — Opred), compared to the observed position angle 
errors ag. We thereby created a map of y 2 values; the x 2 
of a given pixel is a measure of the likelihood that the mea- 
sured polarizations in neightboring pixels have arisen due 
to scattered photons originating in the pixel of interest. 

The mean value of \ 2 was 8.3 x 10 3 ; several minima, 
with x 2 of order 1-10, were present. To minimize the im- 
pact false (noise) minima, we averaged over 6' regions to 
produce our final x 2 map (Figure [p3) . We have highlighted 
the regions of low x 2 by plotting Figure (3b using contours 
at 1,2,3,4, and 5. The smoothing procedure has failed to 
eliminate the worst of the false minima, which are typi- 
cally shallow and spatially compact, or adjoin an edge of 
the image. Broad but shallow minima which correspond 
to sources such as LH 4 or NGC 1755 (see Figure^, Table 
2) have been obscured by the smoothing procedure and do 
not appear in Figure |Jb. 

We correlated the x 2 an d intensity images in order to 
identify the birthplaces of scattered photons. During this 
process, we discovered that the locations of the x 2 minima 
and the centroids of the OB associations were, in general, 
offset from each other. This behavior is expected for light 
scattered from an extended dust layer, inclined to the line 
of sight. For the Nil complex, which is the brightest re- 
gion in the WISP field and lies in a region of low back- 
ground, the offset between the closest x 2 minimum and 
the association centroid was 8.4 j^'® arcminutes in the — x 
direction (roughly west), and 1.8 ±1.8 arcminutes in the 
+y direction (roughly north). The other minima showed 
wide scatter about these values. 

The identifications, coordinates, and a brief description 
of the 9 identified sources of scattered light are given in 



Table 2. We also list the value of x 2 f° r each object; these 
values are indicative of the degree to which the polarization 
patterns near the objects resemble the centro-symmetric 
halo that would be produced by scattering of starlight from 
the central objects from a thin disk of material. 

6 of the 9 regi ons are large OB associa tions or complexes 
as identified by Lucke fc Hodge (1970) . Two of the oth- 
ers are large, very young open clusters; th e last is an H I I 
region associated with a small starcloud (Henize (1956)); 
the ionizing star is lik ely to be the 11th magnitude Bl su- 
pergiant Sk -70° 29 flSanduleak 1970| ). The illuminating 
sources in Table 2 are marked with star symbols in Figure 

Due to our averaging procedure, and because the er- 
rors in position angle do not form a normal distribution, 
our calculated x 2 values cannot be easily interpreted in 
absolute terms. We therefore do not place formal confi- 
dence limits on our detections of scattering halos around 
the illuminating sources. However, the x 2 values listed in 
Table 2 do accurately reflect the relative strengths of the 
detected scattering halos; the symbol sizes in Fig. || have 
been scaled according to these values. 

In sharp contrast to LH 15 and the ass ociations in N 



11, the group labelled B2 (notation from Martin et al. 



1976) corresponds to a broad maximum in x 2 , suggesting 



that the nearby ISM does not support a scattering halo 
around the OB stars within. Similarly, the open cluster 
NGC 1818, which is of a similar age and brightness to the 
scattered light sources NGC 1755 and NGC 1711, shows 
no evidence for a scattering halo. These differences are 
discussed further in jg. 

6. ANALYSIS 
6.1. Correlations with other wavelengths 

The physical interpretation of our UV surface brightness 
and polarization maps is greatly aided by the consideration 
of observations at optical, infrared, and radio wavelengths. 
UV polarization maps probe the structure of the interstel- 
lar medium surrounding young, massive stars. Such ex- 
tremely young stellar populations are often attended by 
the signatures of recent star formation, such as molecu- 
lar gas, dark nebulae and emission from heated dust. The 
dust which both absorbs and scatters the copious UV emis- 
sion from OB associations is distributed within the neutral 
hydrogen clouds of the galactic disk. While massive stars 
are the most luminous, they are dwarfed in number by 
the less massive stars which form along with them as well 
as those which have survived from earlier star formation 
epochs; the combined light of low-mass stars contributes 
to the UV surface brightness of the LMC, and will tend to 
dilute the observed polarization. 

Dust can be traced directly, via maps of the interstellar 
extinction and by mapping the radiation from heated dust 
grains, and indirectly, via maps of the distribution of H I 
and estimates of the gas-to-dust ratio. In the infrared, it 
has been found that the 60 /im and 100 /xm data taken by 
the IRAS satellite correlate well with the H I distribution. 



As discussed by pauvage fc Vigroux (1991 J , the bulk of the 
60-100 /im emission is due to the dust grains in the gen- 
eral ISM. We note a good correlation between the maps of 
Schwcring (1989) and our UV image; the surface bright- 
ness at 60-100 /mi and in our Fig. are of the same order 
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of magnitude, with the IRAS values larger. This confirms 
that Balmer continuum photons provide the bulk of the 
energy that heats the interstellar dust. 



Maps of the interstellar reddening, e.g., Isserstedt & 
are more sensitive to dust above the plane 



(1984) 



Kohl _ 

of the disk than within it, but reveal that high amounts of 
extinction are present in front of N 11, in the area extend- 
ing from LH 17 to LH 20, at the end of the stellar bar, and 
over a large region west of 5 h right ascension and south of 
—67° 30'. There is little apparent correlation between re- 
gions iof_lfijdT_F£B_v_aim^ studied by, 
e.g., ieVaucoulcurs fc Freeman (1972)| . 



The distribution of H I, mapped at low resolu tion by 
Luks & Rohlfs (1992) and at high resolution by Kim e\ 



~al. (1 |998) , correlates closely with the IRAS maps; high- 
column density regions are associated with regions of high 
Eb-v, although the correlation is not perfect. These im- 
ages are plotted for comparison in Figure ^, along with an 
optical image. For the direct UV light of OB associations, 
there is a clear anti-correlation between H I column den- 
sity and UV brightness; extreme examples are the areas 



of the supergiant Ha shells from Meaburn (1980). The 



shell LMC 1 surrounds the OB complex LH 15, and light 
from that association escapes virtually unabsorbed. The 
same is true for the giant shell around N 11; B2 is sur- 
rounded by a network of shells, and LH 12 lies in the wall 
of the supergiant shell LMC 6. The reason for the anti- 
correlation is clear; stellar winds and supernovae clear the 
areas around prominent young associations and clusters. 
However, when regions far from the large associations are 

considered, the picture is no t as clear. 

We examined the maps of Kim et al. (1998), looking for 
correlations between structure in the H I distribution and 
the UV light. Figure |5| shows that instances can be found 
where the diffuse UV light closely traces the H I distri- 
bution. Three of Meaburn's (1980) supergiant Ha shells 
fall within the WISP field of view and are shown in the 
lefthand panel of Figure |^. For example, the wall of LMC 
1, which extends northeast from N 11, shows a high UV 
surface brightness; throughout the image, filaments and 
knots of H I can be identified with UV emission. The lack 
of H I east-southeast of N 11 is seen as well in the UV; the 
same is true for the regions north of NGC 1755, and in 
the extreme southwestern corner of the image. However, 
counterexamples are present as well; the UV spur extend- 
ing from LH 15 to NGC 1818 corresponds to a marked lack 
of H I. 

As shown in Figure ||, the optical emission from the gen- 
eral stellar population provides a crucial piece of the puz- 
zle as well. Due to the smoother distribution of late-type 
stars, the spiral arm structures are much less prominent 
in the optical image. These stars, while faint in the UV, 
nevertheless contribute by their sheer numbers to the UV 
surface brightness. While individual stars are generally 
beneath our detection threshhold, they contribute an un- 
polarized, unresolved background that will tend to lower 
our polarization measures. The stellar density is highest 
in the southeast quadrant of our field, where we see the 
western end of the LMC's bar. The stellar density falls off 
rapidly northward of —68°, an impression confirmed by an 
examination of a starcount map produced from the USNO- 
A astrometric catalog (Levine & Monet, private commu- 
nication). 



It is apparent that we are looking at a field whose prop- 
erties show a strong north-south gradient. In the south, 
the UV light correlates well with general stellar density, 
identifications with features in the H I and IRAS maps are 
not obvious, and potential sources of scattered light, such 
as B2, seem less likely to contribute polarized flux to the 
UV light. In the north, the field stellar density drops, the 
UV light correlates somewhat more strongly with the H I 
surface brightness, and the large, young associations, such 
as N 11 and LH 15, seem to support extended scattering 
halos. 

6.2. Scattering geometry 

Given the general trends with environment explored in 
the previous section, we now examine the effects of the 
placement of illuminators and scatterers on the observed 
polarization map. We wish to uncover the similarities be- 
tween the sources in Table 2; to discover why the B2 com- 
plex and the cluster NGC 1818 seem not to contribute to 
the polarization map; and to assess the possibility that 
stars from beyond the WISP field contribute to the scat- 
tered light. 

The most obvious requirement for the existence of scat- 
tering halos around bright associations is the presence of 
interstellar dust around the associations. The asymme- 
try of the scattering phase function must also be consid- 
ered; there is strong evidence for a forw ard-throwing UV 
phase function for inters tellar dust (e.g., Murthy & Henry 
(1995) , Witt et al. 1992] ). In these cases, back-scattering is 



suppressed, and associations that lie between the observer 
and a scattering layer will produce smaller scattered light 
intensities than associations that lie behind the dust. 
The WISP field of view encompasses two of the spiral 



arms identified by deVaucouleurs fc Freeman (1972 



arm 



A extends from LH 15 down east side of our image to LH 
17 and is defined by a chain of clusters and associations; 
arm B extends from N 11 to the B2 c omplex in the south 
and is less well-defined in our image. [Btecher et al. 1982 



argued that spiral arms rich in OB associations and H II 
regions should produce more direct starlight than scat- 
tered light in the U V. Additionally the lack of dust in the 
LMC's spiral arms (Oestreicher & Schmidt-Kaler (1996)) 



will tend to suppress scattering. 

A common feature of the sources in Table 2 that sup- 
port scattering halos is their placement within or at the 
edge of H I giant or supergiant shells. However, given the 
generally chaotic distribution of H I, this may be a difficult 
situation to avoid. 

One bright UV source that does not support a scattering 
halo is NGC 1818, prominent in the UV spur containin g 
LH 15. NGC 1818 is 20-40 Myr old ( [Westcrlund 1997| ), 
and is brighter than both of the detected clusters NGC 
1755 and NGC 1711 (V 18 is = 9.70, V1755 = 9.85, V m i 
= 10.18; van den Bcrgh 1981). The region surrounding 
NGC 1818 is remarkably quiescent in the H I map, and the 
cluster itself is only m oderately reddene d, with E^-v = 
0.05. The IRAS map ( |Schwering (1989)D also shows little 
thermal dust emission in the NGC 1818 region. The WISP 
UV image shows moderate surface brightness around NGC 
1818, but in the higher resolution image of SCH, it begins 
to be resolved into individual stars. We conclude that 
NGC 1818 is prevented from contributing to the polariza- 
tion map by a lack of nearby dust. 
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Figure 5 
is available in .gif format 
from http://xxx.lanl.gov/ 



Fig. 5. — A multiwavelength view of the western LMC; far Jeft: 6600 A image from a 10 inch refractor (Sandagc 1961); left center: 2150 A 
image from WISP; right center: 60 (im image from IRAS (Schwcring 1989); far right: 21 cm H I image from the ATCA (Kim et al. 1998). 
The locations of supergiant Ha shells (Meaburn 1980) are overplotted on the optical image. 
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The other perplexing association is the large B2 complex 
in the southern WISP field. Consisting of the Lucke-Hodge 
associations LH 1, 2, 5, and 8, B2 is nearly the equal of 
N 11 in the WISP image, and appears to coincide with a 
high column of H I. The field stellar density around B2 
is high, but not as high as around the LH 16-17-20 area; 
these associations are much smaller than B2 but still show 
some evidence of a scattering halo. 

There is clearly material available to scatter the starlight 
of B2: its western half is highly reddened, an unambiguous 
signature of dust. The OB stars of B2 also power several 
H II regions, most notably N 79 and N 83; these regions 
are expected to appear as reflection nebulae in the UV. 
We speculate in passing here that the LH 8 starcloud may 
not be associated with the rest of B2; it is less highly red- 
dened and contains no large H II regions. Moreover, it is 
spatially coincident with the LMC 7 supergiant Ha shell 
of Me abur n (1980), corresponding to a hole in the Kim e\ 
al. (1 398) H I maps. There is another large hole in the H 
I distribution at the western end of the bar (Kim, private 
communication) . 

One possibility is that the B2 complex is above the main 
disk of the LMC, and its scattering halo is therefore sup- 
pressed by the forward-scattering phase funct ion of the 



dust grains. According to the channel maps of Kim et al 
(1998 J, most of the H I around B2 lies at a heliocentric 



velocity of 241-274 km s _1 . Wayte (1990) has stated that 



interstellar absorption lines due to LMC gas are present in 
stars in this region at velocities of «220 km s _1 , but not 
at higher velocities, implying that the illuminating stars 
lie above the main LMC disk. 

A further possibility is that the large angular extent of 
B2 exceeded the limits of the window we used to analyze 
X 2 in the polarization map. The 42 arcminute-square win- 
dow only just contained all of B2, and the point-source 
approximation we used to test against the distribution of 
position angles was almost certainly invalid for this case. 
In addition, most of the window was probably swamped 
with direct stellar light. We repeated the x 2 analysis with 
54 and 63 arcminute windows, but the results for B2 re- 
mained unchanged. The persistent non-detection leads us 
to suggest that B2 does in fact lie above the main LMC 
disk, but this is highly uncertain. 

The results of our x 2 test suggest that scattering of light 
from point sources by a smooth disk cannot be the only 
source of polarization in the WISP field. In this case, it 
is possible that light is being scattered from outside the 
WISP field of view, perhaps from many kiloparsecs away, 
into our line of sight. Prime candidates for illuminating 
sources would be the OB associations in Shapley Constel- 
lation IV (50' east of LH 19), and the associations in the 
bar, extending southeast of L H 20. Further afield, t he 30 
Dor association, postulated by Luks fc Rohlfs (1992)| to lie 
as much as 400 pc above the LMC midplane, could be a 
contributor to scattered light in the WISP field. 

To be able to scatter light from associations across the 
LMC, dust clouds must be raised above the plane of the 
disk, which is optically thick. Among L* spirals, galax- 
ies with comparable rates of star formation to the LMC 
have often had large amounts of dust ejec ted to distances 
a kiloparsec or more above the plane (e.g., Howk & Savage 
1997|] Given the propensity of spiral galaxies for ejecting 
dust into their halos, and the previously suggested depth 



structure of the LMC, it is quite likely that the LMC pos- 
sesses extraplanar dust; some of this dust may scatter the 
light of distant OB associations. 

7. DISCUSSION 

7.1. To What Extent is the Western LMC a 
Kiloparsec-Scale Reflection Nebula? 

Several authors have suggested that UV-bright, gas rich 
galaxies possess extended halos of scattered UV light, in 
effect creat i ng galaxy-sized r eflection nebulae (e.g., Stechci 
et al. 198S, |Neff et al. 1994 ). Characteristics of scattered 
light halos should include somewhat bluer spectral energy 
distribution than the overall radiation field of the galaxy, 
a correlation with the H I distribution, and linear polar- 
ization. Does the western side of the LMC disk meet any 
of these criteria? 

The WISP field is richer in young stars than the far 
outer regions of spirals studied by Stecher et al (1982) 
and Neff et al. 1994. Thus any putative scattered light 
contribution to its UV surface brightness will be diluted 
by direct, stellar light. In the southern half of the WISP 
field, the UV surface brightness correlates closely with op- 
tical images and with star count maps, suggesting that 
direct light accounts for most of the ultraviolet radiation. 
However, northward of « —68°, the stellar density falls off 
rapidly; Figure || shows that the ultraviolet light falls off 
more slowly. Throughout the image, there are areas of dif- 
fuse UV emission that begin to be resolved in the images 
of SCH, but in the northern WISP field it is easy to see 
correlations between diffuse UV emission and the presence 
of H I emission. 

The degree to which scattered light contributes to the 
total UV surface brightness of the LMC can be assessed 
by examining the mean level of polarization. Our error- 
weighted mean P = 12.6% ±2.3%; in galactic reflection 
nebulae, the polarization at optical wavel engths is typi- 
cally ~20-40% (com parable to p ma x', e.g., IWatkin, Gled- 



hill & Scarrott 1991). We assume that the pixels in which 



wc find P = 30-40% are analogous to reflection nebulae; 
light that has been scattered nearly in the plane of the sky 
is the primary contributor to their flux. However, such 
high values are only found in areas of low flux, and thus 
high error, in our map. 

Because the dust scale height of the LMC disk is non- 
negligible, and because of the inclination of the disk to 
the plane of the sky, typical polarization levels should fall 
below Pmax- Scarrott et al. (1990) showed the variation of 
P with scattering angle for several dust grain size distri- 
butions. Polarization levels of 10-15% are to be expected 
for scattering angles B « 50°-60° (and for 9 « 125°-135° 
if backscattering is important). These angles are close to 
the complement of the LMC's angle of inclination to the 
sky (28°-45°), and would be reasonable values for a mean 
angle of scattering. 

As noted in most of our measured polarization values 
lie between 5% and 20%; this implies that in the average 
region of the WISP field, at least 20%, and possibly >90% 
of the detected light has been scattered. The statistical 
significance of many of the polarization vectors is low; we 
seem able to reliably detect polarizations where the con- 
tribution of scattered light is ;>50%. 

Jui'a (1980)' presented a method to predict the intensity 
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of scattered light from a dust cloud at a certain angular 
distance from an illuminating source. We adapted their 
method to the case of forward-scattering grains; this model 
predicts that scattered light surface brightnesses near com- 
plexes such as Nil and B2 should lie in the range 0.5 - 
1.5 x 1(T 8 erg s" 1 cm" 2 A^ 1 Sr" 1 , within 30' of the il- 
luminating source s. This level should be taken as a strict 
lower limit, since Jura (1980) considered dust clouds in 
the plane of the sky. Accordingly, the WISP data show 
surface brightnesses some 4-20 times higher than the Jura| 
(1980fl formula predicts. While some of this discrepancy is 
due to the integrated contribution of scatterers along the 
line of sight (in front of or behind the plane of the sky) , it 
does suggest that scattered light is not the only contribu- 
tor to the total surface brightness. The largest difference 
between the predicted intensity due to pure scattered light 
and the observed surface brightness occurs for the south- 
ern part of the WISP field, where the field stellar density 
is high. 

We conclude from the mean level of polarization and 
the distributions of UV surface brightness around OB as- 
sociations that scattered light contributes at least 20% to 
the total UV surface brightness outside the associations, 
and possibly more than 50%. The scattered fraction is 
expected to be higher in the north, where stellar densities 
are lower and we more confidently detect centro-symmetric 
polarization patterns around the largest associations. In 
this respect, the northwestern disk of the LMC can be re- 
garded as a very faint UV reflection nebula. The referee 
has suggested that similar observations in another band- 
pass could be a test of this hypothesis. Because of the dif- 
fering spectral energy distributions of the reflected light of 
early-type stars and the unresolved sea of late-type stars, 
an ultraviolet color map could test the idea that the scat- 
tered light fraction is higher in the north. Observations at 
1500 A could be obtained, although the ability to obtain 
sufficient signal-to-noise is in question. 

7.2. What are the Scattering Geometry and Phase 
Function? 

The geometry of the scattering material is of crucial im- 
portance for determining the distribution of polarizations. 
We have seen that the B2 complex shows no evidence for a 
scattering halo, and suggested a possible geometric expla- 
nation. Another issue is the offset of centro-symmetric po- 
larization patterns from the central OB associations, noted 
in §|[ A possible explanation for this offset is the inclina- 
tion of the LMC disk to the line of sight. The LMC is 
known to be tilted between 28° and 45° from the plane of 
the sky, along a roughly north-so uth axis, with the east- 
ern (30 Dor) side closer to Earth ( |Westerlund 1997| ). The 
effects of varying the inclination of a scattering disk ha ve 



been investigated by Cole, Wood, & Nordsieck (1999b) 



this context, and this explanation does seem plausible. 

Another scattering issue that requires modelling is the 
asymmetry in phase function: do the dust grains scatter 
isotropically, or are they strongly forward-throwing? If 
isotropic, the UV surface brightness should fall off more 
slowly than for the forward-throwing case; the mean po- 
larization level may also change as the preferred paths for 



photon escape from the disk are changed. For isotropically 
scattering grains, the geometric explanation for B2's lack 
of scattering halo is weakened. 

8. SUMMARY 

We have obtained the first wide-field ultraviolet polari- 
metric images in astronomy, using the rocket-borne Wide- 
Field Imaging Survey Polarimeter (WISP). We find the 
diffuse UV background light is polarized at levels from 4- 
40%, the lower limit set by our sensitivity threshold. 

We photometrically c alibrated our data using th e data 
of [Smith et al. (1987) for the OB associations of Luckc 
fc Hodge (1970)| ; our photometric accuracy is probably ± 
«20%. The UV surface brightness (Fig. g) is of the same 
order of magnitude as the IRAS 60/im and 100/im mea- 
sures, indicating that Balmer continuum photons are the 
primary heaters of interstellar dust in the LMC. We find 
a minimum UV surface brightness of 5.6 ±3.1 x 10~ 8 erg 
s _1 cm -2 A -1 Sr -1 (23.6 ±0.5 mag arcsec -2 ). 

A polarization map (Fig. |^a) was created and analyzed 
in order to study the diffuse galactic light of the LMC. The 
mean level of polarization was P = 4.4% ± 7.8%; however, 
we excluded regions where P < 2.5cr from the polarization 
map; the mean of the statistically significant polarization 
vectors was P = 12.6% ±2.3%. This suggests a mean con- 
tribution of scattered light to the total of >20%. We used 
a modified % 2 test to identify the sources of scattered light 
within the WISP field, and found 9 bright regions that 
might be responsible for scattered light. These are listed 
in Table 2 and Figure |[ 

The fall-off of surface brightness with radial distance 
from the OB associations was examined and suggests that 
the scattered light is not the sole contributor to the ob- 
served surface brightness. There is a north-south gra- 
dient of properties across the field, due to the falloff in 
stellar density with distance from the bar. The north- 
ern half of the image contains the strongest evidence of 
centro-symmetric polarization patterns and large amounts 
of scattered light contribute to the UV surface brightness 
above -68°. 

We find that the B2 complex shows no evidence for a 
contribution to the scattered light around it; it does not 
support a scattering halo. This may be an effect of the 
high field stellar density in the neighborhood and unfavor- 
able geometric positioning of stars relative to dust. Al- 
ternatively, it may be the result of bias introduced in our 
analysis of the polarization map. 

Many of the questions raised by these data can only be 
answered at present through modelling efforts. Some of 
the issues of sca ttering geometry and phase functio n will 
be addressed by Cole, Wood, & Nordsieck (1999b) in an 
accompanying paper. 
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